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ABSTRACT 

We report the results of a multi-wavelength study of the nearby galaxy group, Abell 3581 
(z = 0.0218). This system hosts the most luminous cool core of any nearby group and exhibits 
active radio mode feedback from the super-massive black hole in its brightest group galaxy, 
IC 4374. The brightest galaxy has suffered multiple active galactic nuclei outbursts, blowing 
bubbles into the surrounding hot gas, which have resulted in the uplift of cool and cold gas into 
the surrounding hot intragroup medium. High velocities, indicative of an outflow, are observed 
close to the nucleus and coincident with the radio jet. Thin dusty filaments accompany the 
uplifted, ionised gas. No extended star formation is observed, however, a young cluster is 
detected just north of the nucleus. The direction of rise of the bubbles has changed between 
outbursts. This directional change is likely due to sloshing motions of the intragroup medium. 
These sloshing motions also appear to be actively stripping the X-ray cool core, as indicated 
by a spiraling cold front of high metallicity, low temperature, low entropy gas. 

Key words: Galaxies: groups: individual: Abell 3581 - galaxies: clusters: intracluster 
medium - galaxies: active - galaxies: ISM 



1 INTRODUCTION 

The last decade has seen significant progress into the cooling flow 
phenomenon in 'cool core' clusters of galaxies. These clusters are 
so named as they exhibit sharply peaked X-ray surface brightness 
profiles resulting from their cool, high density X-ray cores. The 

* Based on observations obtained at the Southern Astrophysical Research 
(SOAR) telescope, which is a joint project of the Ministerio da Ciencia, 
Tecnologia, e Inovacao (MCTI) da Republica Federativa do Brasil, the U.S. 
National Optical Astronomy Observatory (NOAO), the University of North 
Carolina at Chapel Hill (UNC), and Michigan State University (MSU). 
f E-mail: rcanning@stanford.edu 



cooling time of the gas in these cores is significantly shorter than 
the Hubble time. 

If no sources of heat are present, the observed high central 
X-ray luminosities would imply mass cooling rates as high as hun- 
dreds or thousands of solar masses a year (e.g. Fabian 1994 and 
references therein). However, with the advent of the XMM-Newton 
space telescope, high resolution spectra exposed a lack of strong 
X-ray emission below 1 keV in cool cores clusters. The implica- 
tion of these observations was a reduction in the inferred net X-ray 
cooling rates and an effective quenching of any resulting star for- 
mation (for reviews see Peterson & Fabian 2006 and McNamara & 
Nulsen 2007). The most probable source of heating to quench the 
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Figure 1. Left: Three colour image of Abell 3581 formed from an adaptively smoothed Chandra X-ray image with the average at each radius removed (blue), 
archival VLA 1.4 GHz radio emission (green) and narrowband Ha+[N n] image (red). The radio bubbles push aside the hot X-ray gas forming the inner 
two cavities. Another 'ghost' X-ray cavity is observed to the south-west devoid of radio emission. Filaments of ionised gas encase the young bubbles and are 
entrained beneath the older cavity. Right: Colour composite of the central group galaxy, IC 4374, in the three HST filters F300X (blue), F555W (green) and 
F814W (red). The region shown to the right is highlighted on the left hand figure by the white rectangle which has dimensions 24.2 arcsec (10.5 kpc) across. 
Dust filaments are observed aligned with the Hcn+[N n] and a young star cluster is observed in the dust to the north of the nucleus. 



cooling of the intracluster medium (ICM) is provided by outbursts 
from the super-massive black hole (SMBH) lurking at the centre 
of the central, dominant, brightest cluster galaxies (BCGs) (for a 
review of the evidence for BH feedback see Fabian 2012). 

The prevalence of Ha emission in BCGs has long been known 
(e.g. Lynds 1970; Cowie et al. 1983; Hu et al. 1983; Johnstone et al. 
1987; Heckman et al. 1989; Crawford et al. 1999; Donahue et al. 
2000; Fabian et al. 2001; Conselice et al. 2001) and has been shown 
to have diverse and dramatic morphologies. These ionised gas neb- 
ulosities surround about a third of all BCGs but are found distin- 
guishably in cool core clusters. The luminosity of the Ha bright 
nebulosities ranges from 10 39 — 10 43 erg s _1 and can extend over 
70 kpc from the core of the BCG (Crawford & Fabian 1992; Craw- 
ford et al. 1999). The emission is well correlated with both the X- 
ray luminosity and with emission lines from molecular gas. Devel- 
opments in optical astronomy, particularly with the progression of 
integral field spectroscopy (IFS), have allowed the spatio-kinematic 
properties of the gaseous nebulae surrounding BCGs to be better 
probed. For example, the spatial, kinematical and excitation prop- 
erties of the ionised and molecular gas has now been mapped in 
small samples of BCGs (e.g. Wilman et al. 2006; Hatch et al. 2007; 
Wilman et al. 2009; Edwards et al. 2009; Farage et al. 2010; Oonk 
et al. 2010; Canning et al. 2011). 

In addition to ionised gas nebulosities, many cool core BCGs 
contain significant masses of cold molecular gas and dust, both 
in the centres of the galaxies and farther out. This is often coin- 
cident with the extended ionised gas filaments in the galaxy out- 
skirts (e.g. Jaffe & Bremer 1997; Donahue et al. 2000; Edge 2001; 
Salome & Combes 2003; Hatch et al. 2005; Salome et al. 2006; 
Johnstone et al. 2007; Oonk et al. 2010; Edge et al. 2010; Mittal 
et al. 2011; Werner et al. 2012). What the origin of these cool and 
cold gas phases are and how the complex heating and cooling bal- 
ance in these systems is maintained requires observations of mul- 



tiphase material; from synchrotron emitting relativistic plasmas to 
cold molecular gas and dust. 

Groups and poor clusters have shallower potential wells com- 
pared with rich galaxy clusters. An important implication of this is 
that the mechanical output from an Active Galactic Nucleus (AGN) 
outburst is of the same order as the group binding energy. The ef- 
fects of AGN feedback in groups may well be different to that ob- 
served on the cluster scale. The majority of galaxies in the Universe 
are situated in group environments. Therefore, investigating the na- 
ture of AGN feedback in lower mass systems is paramount to our 
understanding of galaxy evolution. 

The nearby galaxy group, Abell 3581 (Fig. 1, z = 0.0218, 
distance d = 93.6 Mpc), presents an excellent case study for the 
investigation of the thermal regulation of groups of galaxies. It has 
a system temperature of ~ 1.7 keV (Sun et al. 2009) and hosts the 
most luminous X-ray cool core among galaxy groups with z ^0.05. 
The central AGN, PKS 1404-267, has the highest radio power of 
any nearby group with a large cool core (e.g. Sun 2009). This cool 
core co-exists with one of the most spectacular and luminous Ha 
nebulae observed in galaxy groups, which surrounds the brightest 
group galaxy (BGG) IC 4374 (Johnstone et al. 1987; Danziger & 
Focardi 1988). CO(2 — 1) emission is also detected in the central 
galaxy (Edge, private communication). The previous Chandra ob- 
servations reported in Johnstone et al. (2005) have shown that two 
X-ray cavities exist in the centre of the system, where the AGN jets 
are currently displacing the hot interstellar medium (ISM). Sanders 
et al. (2010) observed Fe XVII emission lines in the core of Abell 
3581, indicating the presence of low temperature, 0.5 keV X-ray 
emitting gas. The authors find evidence that this X-ray coolest de- 
tected gas phase has a low volume filling fraction (3 per cent) so 
the cool blobs exist interspersed in the hotter surrounding medium. 
Such low volume filling factor cool X-ray gas has been seen in 
other objects (e.g. Sanders & Fabian 2007; Sanders et al. 2009, 
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Table 1. Summary of observations. 



Telescope 


Filter 


Exposure 


Date 


ID 


PI 


SOAR/SOI 
HST/WFC3/UVIS 

VLT/VIMOS 


CTIO 6693/76 
CTIO 6520/76 

F300X 

F555W 

F814W 
HRO grism, 0.67" fibers 


3x600 s 
3x360 s 
5427 s 
1335 s 
1044 s 
4500 s 


2008 July 8 
2008 July 8 
2011 February 2nd 
2011 February 2nd 
2011 February 2nd 
2009 March 30th 


12373 


Sun/Donahue 
Sun 

Fabian 


EVLA 
GMRT 


L band, BnA configuration 
610 MHz 
235 MHz 


4hr 
2.4 hr 
2.4 hr 


2012 September 8th 
2010 May 
2010 May 


SB0410 
18-009 


Sun 
Sun 


Chandra/ACIS-S 




85 ks 


2011 January 3rd 


12884 


Sun 



2010; Werner et al. 2010, 2012). Sanders et al. (2010) find no ev- 
idence for O VII emission lines in Abell 3581 implying very little 
of the gas can be cooling radiatively below 0.5 keV. Either non- 
radiative cooling is occurring or a heating source is needed which 
targets this coolest X-ray gas without overheating the surroundings. 
Thus, Abell 3581 is an ideal system in which to study AGN heating 
and the multi-phase gas in group cool cores. 

In this paper we present deep multi-wavelength imaging 
and spectroscopy of the central galaxy in Abell 3581 and the 
surrounding ICM (details of the observations are given in Table. 1). 
Throughout this paper we assume the standard ACDM cosmology 
where H = 71 km s 1 , Q m = 0.27 and Qa = 0.73. For this 
cosmology and at the redshift of IC 4374 (z — 0.0218), an angular 
size of 1" corresponds to a distance of 0.435 kpc. 

2 OPTICAL DATA 
2.1 SOAR 

Narrow-band optical imaging was performed with the 4.1m South- 
ern Astrophysical Research telescope (SOAR) telescope on July 
8, 2008 (UT) using the SOAR Optical Imager (SOI). The night 
was photometric with a seeing of around 1.0". The CTIO 6693/76 
narrow-band filter was used for the Ha+[N II] lines, and the CTIO 
6520/76 filter for the continuum. Three 600 sec exposures were 
taken with the 6693/76 filter and three 360 sec exposures with the 
6520/76 filter. Each image was reduced using the standard proce- 
dures in the IRAF MSCRED package and the spectroscopic stan- 
dard EG 274. The pixels were binned 2x2, for a scale of 0.154 arc- 
sec per pixel. More details on the SOI data reduction and the con- 
tinuum subtraction can be found in Sun et al. (2007). The net (con- 
tinuum subtracted) Ha image has been published in Sun (2012). 



2.2 HST 

Observations were made using HST/WFC3/UVIS on 2011 Febru- 
ary 2 in the F300X, F555W and F814W filters. The total expo- 
sure times in each filter were 5427, 1335 and 1044 s respectively. 
Throughput of the F555W filter at the redshifted wavelengths of 
the strong emission lines, Ha+[N II] ^6700 A, is only of order a 
few percent. Standard procedures were used to bias- subtract and 
flat-field the data frames and to drizzle the science frames together, 
these procedures use the software DrizzlePac and are described in 
the HST WFC3 data handbook. 



The optical images are shown in Figs. 2, 3 and 4. No cor- 
rections for reddening have been made to the images, although a 
relative intrinsic E(B — V) image of the galaxy is shown in the 
right-hand panel of Fig. 4. This reddening map is formed by first 
correcting the broad band HST imaging for Galactic reddening in 
the direction of Abell 3581. We take the Galactic reddening to be 
E(B-V) = 0.060 ± 0.001 using the reddening maps of Schlegel 
et al. (1998) and assuming the Galactic extinction law of Cardelli 
et al. (1989). We then converted the HST/WFC3/UVIS broad band 
images into AB magnitude units and subtracted the F814W broad 
band HST image from the F555W band. We then subtracted the 
value at a reference point, within the galaxy, but where there was 
no evidence for dust lanes in the single F555W or F300X frames. 
The point chosen as our reference is given by the coordinates RA 
14^07 m 29 s .768, Dec. -27° 01' 05.20" (J2000). This then pro- 
vides an E(F555W — F814W) relative intrinsic reddening map. 

To convert the map to E(B — V) units requires the assump- 
tion of an attenuation curve. We assume the Cardelli et al. (1989) 
extinction law and calculate A\/Ay from the HST filters. We use 
the pivot wavelengths of the F555W and F814W filters given in the 
HST WFC3 data handbook as the effective wavelength of the fil- 
ters. Due to the broad nature of the filters and the low redshift of 
our object, no K-correction is made. 



2.3 Integral field spectroscopy 

Observations were made on 2009 March 30 using the Visible Mul- 
tiObject Spectrograph (VIMOS) on the VLT in La Paranal, Chile 
(see Le Fevre et al. 2003 and Zanichelli et al. 2005 for a description 
of the VIMOS IFU and a discussion of data reduction techniques). 
We obtained High Resolution Orange (HRO, R = 2650) data using 
the larger 0.67" fibers giving a field-of-view of 27"x27" with 1600 
fibers. The wavelength coverage across the whole field-of-view in 
the HRO grism is from 5250 - 7400 A which cover the redshifted 
spectral lines of [N l]A5199 to [S ll]A6731, including the strong 
spectral lines of [O l]A6300, Ha and [N ll]A6548,6583. However, 
due to the centring of the pseudo- slits in each quadrant the spec- 
tral coverage varies from quadrant-to-quadrant. The result is that 
the eastern side of our field-of-view has additional coverage of the 
redshifted spectral line of [O lll]A5007. The total exposure time of 
these observations was 4500 seconds. 

The data were reduced by the VIPGI 1 pipeline (see Scodeg- 

1 VIMOS Interactive Pipeline Graphical Interface (VIPGI), obtained from 
http://cosmos.iasf-milano.inaf.it/pandora/. 
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gio et al. 2005 for a description of VIPGI). The data cubes were 
combined, transmission corrected, sky absorption and emission 
corrected and corrected for Galactic extinction (E(B — V) — 
0.060 ± 0.001) using a set of IDL routines. The 2 telluric ab- 
sorption feature at ~ 6800 - 6900 A falls over the [S II] doublet 
and as such warranted special attention. We use observations of two 
standard stars taken just before and after the science observations 
to correct for this feature. 

Due to the nature of the filamentary system we choose to bin 
the data based on the surface brightness of [N ll]A6583 emission, 
the strongest line observed, using the contour binning technique of 
Sanders (2006) and binning to a signal-to-noise of 10. 

We then fit the emission lines in IDL using MPFIT (More 1978; 
Markwardt 2009). The strong emission lines of [N ll]A6548,6583 
and Ha are fit simultaneously across the field-of-view. The redshift 
and velocity dispersion are constrained to be the same for these 
emission lines and the integrated flux of the [N II] doublet is tied 
in the ratio 3:1, the scaling being dictated by the atomic parame- 
ters (Osterbrock & Ferland 2006). The [S II] and [O I] lines are 
fit separately (the [O I] doublet flux is also tied by atomic parame- 
ters); however, the morphology and velocity structure of all emis- 
sion lines are similar. We fit gaussians to the spectra with 1 and 2 
velocity components and both with and without a broad Ha com- 
ponent. Using MPFTEST in IDL we determine that a one component 
velocity structure is sufficient to fit the majority of our data. How- 
ever, we require a broad Ha component within the inner 2 square 
arc seconds radius of the galaxy. 

In the images from IFU data, only pixels where the emission 
lines were detected at a 4 a or greater level are shown. The instru- 
mental spectral resolution is determined from a fit to strong sky 
emission lines and is found to be 98 krns" 1 at full width half 
maximum (FWHM). 



3 X-RAY DATA 

Chandra X-ray data of Abell 3581, with an exposure time of 85 ks, 
were taken on 2011 January 3rd (ObsID 12884). The observa- 
tions were conducted with the Advanced CCD Imaging Spectrom- 
eter (ACIS). The data was reprocessed using the CI AO (Fruscione 
et al. 2006) ACIS -PROCESS .EVENTS tool, applying the VFAINT 
processing to lower the non X-ray background. To look for contam- 
inating flares, we examined the lightcurve of ACIS CCD 5 between 
2.5 and 7 keV in 200s bins. A sigma clipping algorithm was used 
to search for flared periods, of which there were none. Standard 
blank- sky background datasets were used, re-projected to match 
the observation and with the exposure time adjusted to match the 
background-dominated 9 to 12 keV observational count rate. To 
make the spectral maps, we used the Contour Binning algorithm 
(Sanders 2006) to create regions with a signal to noise ratio of 
20 (around 400 counts). Spectra and background spectra were ex- 
tracted from these bins, and response and ancillary response regions 
were created, weighting response regions by the number of counts 
in the 0.5 to 7 keV band. We fitted the spectra in Xspec between 0.5 
and 7 keV, minimising the modified C statistic. We used the APEC 
thermal model to fit the spectrum (Smith et al. 2001), with Galactic 
absorption (4.36 x 10 20 cm -2 ) applied. In the fits the temperature, 
metallicity and normalisation were allowed to vary. 

In order to examine the surface brightness edges in the images 
we fit and subtracted a double-/^ model, using the Sherpa fitting 
package, to an azimuthally averaged, background subtracted sur- 
face brightness profile, centred on and extracted between a radius of 



2.5-120 arcsec. The surface brightness profile is from the exposure- 
corrected image with point sources removed manually. The image 
is in the energy range 0.5-7 keV. Despite the asymmetry of the clus- 
ter, the model fits the data reasonably well. However, deviations 
from the fit are seen at ^50 arcsec, ~80 arcsec, and in the central 
20 arcsec. These deviations correspond to surface brightness edges 
which will be examined further in Section 5.2. Double-/^ fits to 
individual sectors and empirical azimuthally averaged profiles are 
also fit and removed from the surface brightness images to high- 
light substructure. 

For the temperature and density profiles across the 'bubble' 
and 'cold front' (see Section 5.2) we extracted spectra from the 
regions of interest, and fitted the spectrum within Xspec (version 
12.8.0a 2 ) with an absorbed MEKAL model. The regions for the 
spectral analysis were chosen such that there were at least 3000 
counts in each bin. Only data in the 0.5-7.0 keV range were used. 
C statistics were used for the fit and spectra were deprojected us- 
ing the routine developed in Sanders & Fabian (2007) and Russell 
et al. (2008). The Galactic absorption in the direction of Abell 3581 
was taken as 4.36 x 10 20 cm -2 , using the value of Kalberla et al. 
(2005) and was frozen in the fits. The solar photo spheric abundance 
table by Anders & Grevesse (1989) was used in the spectral fits. Er- 
ror bars on the temperature, density and pressure plots are shown at 
the 90 per cent level. 

4 RADIO DATA 
4.1 EVLA 

High frequency radio observations of Abell 3581 were obtained 
with the NRAO 3 EVLA (Perley et al. 2011) at L band (1-2 GHz). 
Due to the low declination of the target, the observations were un- 
dertaken in the hybrid BnA configuration on 2012 Sept. 08 (project 
ID: SB 04 10). The extended north arm of the hybrid configuration 
allows for a more circular beam for the target. The observations 
were a total of 4 hours and were taken with two 512 MHz wide 
basebands centered at 1314 and 1686 MHz. We configurated the 
WIDAR correlator so that each baseband had 8 subbands with 64 
channels per subband for a frequency resolution of 1 MHz. We used 
3C286 as our absolute flux and bandpass calibrator and J1351-1449 
as our phase calibrator. 

The data were flagged and calibrated using the OBIT data re- 
duction package (Cotton 2008). We Hanning smoothed the data to 
reduce the effects of Gibbs ringing around strong radio frequency 
interference (RFI) which is common at L band. We then flagged 
the data in the time and frequency domains to remove the major- 
ity of the significant RFI. Additional time-domain flagging on the 
bright calibrators was done using the average calibrator amplitudes. 
We used all calibrators to calculate the parallel hand group delays 
and applied the results to all sources. We then followed standard 
techniques for bandpass calibration as well as amplitude and phase 
calibration. Following calibration, we split off the target field and 



2 We used the AtomDB 2.0 release which includes changes to the iron L- 
shell data affecting kT < 2 keV plasma which has the effect of increasing 
the temperature (10-20 per cent) while decreasing the abundances (~20 per 
cent) when compared with AtomDB 1.3.1. 

3 The National Radio Astronomy Observatory is a facility of the National 
Science Foundation operated under cooperative agreement by Associated 
Universities, Inc. 
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Figure 2. Left: Net narrow band Hen emission image of the BGG in Abell 3581 after removal of continuum emission with the CTIO 6520/76 filter. Right: 
The unsharp masked Hen emission image. The initial image has been smoothed by a 2D gaussian with a = 1 arcsec, subtracted from the original and binned 
by 4 pixels (0.6") in each direction. 




Figure 3. Left: HST F555W image with a smooth elliptical isophotal model subtracted. The dark feature piercing north-south through the nucleus and turning 
towards the south-east is absorption from a prominent dust lane. Right: F300X image of the same region. The red circle shows a young star cluster not seen 
in the visible (F555W and F814W) images. The star cluster coincides with the northern dust lane. 




Figure 4. Left: The central region of the net Hen image overlaid with contours of dust absorption. The dust closely coincides with the inner Hen filaments. 
Right: The relative dust attenuation in the brightest cluster galaxy in units of E(B — V) magnitudes after correction for Galactic reddening. Contours of 
L-band EVLA radio emission, which fill the inner X-ray cavities, are overlaid. The angular scale is the same as the left hand panel. 
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undertook additional flagging before imaging the data with MFI IM- 
AGE which allows us to deal with the wideband data. We applied 
one round of phase only self-calibration on the target using our ini- 
tial wide-field wideband image. The final image of Abell 3581 has 
a noise of a=42 /iJy/bm with a resolution of 4.65 x 3.73 arcsec. 



4.2 GMRT 

Abell 3581 was observed with the GMRT at 235 MHz and 610 
MHz (project 18-009) in May 2010, for approximately 3 hours, in- 
cluding calibration overheads. The observations were made using 
the dual-frequency mode. The 610 MHz data were collected using 
the upper and lower side bands (USB and LSB) with an observing 
bandwidth of 16 MHz each. Only the USB was used at 235 MHz, 
with a bandwidth of 8 MHz. We chose the default spectral-line ob- 
serving mode, with 128 channels for each band at 610 MHz (64 
channels at 240 MHz) and a spectral resolution of 125 kHz/channel. 
The datasets were calibrated and reduced using the NRAO Astro- 
nomical Image Processing System package (AIPS) as described in 
Giacintucci et al. (2011). Self-calibration was applied to the data 
to reduce residual phase variations and improve the quality of the 
final images. Due to the large field of view of the GMRT at low 
frequencies, we used the wide-field imaging technique at each step 
of the phase self-calibration process, to account for the non-planar 
nature of the sky. The final images were produced using the multi- 
scale CLEAN implemented in the AIPS task IMAGR, which results 
in better imaging of extended sources compared to the traditional 
CLEAN (e.g., Wakker & Schwarz 1988; Cornwell et al. 1993). The 
full resolution of our final images is 8" x 4" at 610 MHz and 
16" x 11" at 235 MHz. The rms noise level (la) in these images is 
0.1 mJy beam -1 and 0.9 mJy beam -1 , respectively. Residual am- 
plitude errors are within 8 per cent at 240 MHz and 5 per cent at 
610 MHz. 

5 RESULTS 
5.1 Optical 

5.1.1 Morphology 

The BGG in Abell 3581 has a smooth, relaxed appearance in broad 
band optical images. However, the narrow band Hoh-[N II] image 
shows striking filaments of ionised gas. These filaments connect to 
the nucleus in the easterly and westerly directions, parallel to the 
current radio axis, then curve sharply towards the north and south 
respectively (see Fig. 2). The filament system is complex, formed 
of at least a few strands in each direction and is not one continuous 
structure. Filaments of this nature have now been observed in many 
massive galaxies dominating the centres of X-ray bright cool core 
groups and clusters. Our VLT VIMOS IFS data show the emission 
lines of Ha, [N II], [S II] and [O I] in the inner filament all share 
the same morphology. 

A dust lane stretching north to south through the nucleus is 
clearly visible in all optical bands and in the F300X image. Fig. 3 
shows the F555W image with the smooth elliptical component sub- 
tracted, highlighting the myriad dust lanes. The nuclear component 
of the dust lane is ^5 arcsec in length and extends out radially in 
two plume-like structures; one to the north and one to the south. 
In the south, beyond a few arcsec, the dust lane abruptly changes 
direction and traces the ionised gas filament to the east (see Fig. 4). 
The dust extends at least 10 arcsec along the filament; whether the 



dust stops here or continues all the way along the filament is un- 
known. A similar extension to the east is seen in the dust emission 
north of the nucleus. Contours of the dust emission are overlaid on 
the net Ha+[N II] image in the left hand panel of Fig. 4. A young 
star cluster is found coincident with the northern dust lane but no 
other obvious young stellar structures are observed in the filaments. 

The right-hand panel of Fig. 4 shows an E{B — V) map 
of the relative intrinsic extinction in the central galaxy due to the 
dust lanes; the reddening due to absorption by dust is significant in 
the central galaxy. This is overlaid with contours of the EVLA L- 
band radio emission from the inner X-ray cavities (Johnstone et al. 
2005). The dust lane in the south and extending towards the east 
of the nucleus also coincides also with the outer edge of the radio 
bubble. 

5.1.2 Nebula kinematics 

Our VIMOS IFS data (Fig. 5) provides evidence that a broad Ha 
emission line is necessary in the central 22 pixels, a region ~2 arc- 
sec in radius (~3 pixels or ~1 kpc radius) and coinciding with the 
largest surface brightness of high frequency radio emission. Our 
PSF is ~2. 1 pix FWHM so the emission should be resolved. The 
broad line region of the AGN is expected to be significantly smaller 
than this, being typically <0.1 pc or for the most luminous quasars 
~0.25 pc. This broad component has a best-fit line-of-sight velocity 
redshifted ~ 150 kms -1 with respect the the narrow Ha emission 
line component within the same region. The breadth of the line is 
~2900 kms -1 and the Ha flux of the broad component in the 
inner 2 arcsec radius is 2.2 1 ±0.03 x 10 -14 erg cm -2 s -1 . 

The line-of-sight velocity and velocity width of the narrow 
component of the ionised gas nebulae is shown in Fig. 6. The units 
of both figures are km s -1 and the line-of-sight velocity zero point 
is taken to be the nucleus (14^07 m 29 s .8, -27°01'06.4") as de- 
termined from the L-band EVLA radio images. Overlaid are the 
contours from the narrowband Ha image. A smooth line-of-sight 
velocity gradient is observed through the nucleus, perhaps indica- 
tive of rotation. We note that no evidence for rotation of the stellar 
component, observed through the NaD absorption line, is seen. 

The rotation observed in the galaxy core extends smoothly 
along the filament, consistent with a common origin of the fila- 
ment gas and the rotating gas in the core. The velocities are low 
in the filaments with line-of-sight velocities of ~50kms -1 and 
velocity dispersions of ~ 150 kms -1 . The kinematics are there- 
fore very similar to those observed in filament systems surrounding 
other BCGs (e.g. Hatch et al. 2005; Canning et al. 2011). 

5.1.3 Emission line fluxes and ratios 

Our HRO data span the redshifted spectral range of [N l]A5199 
to [S ll]A6731 across the whole field-of-view and the redshifted 
spectral range of [O lll]A5007 to [S ll]A6731 in the south-east and 
north-east quadrants. These data are fit within IDL as described 
in Section 2.3. At the resolution of our observations, all emission 
lines detected in our data have the same spatial morphology and, 
aside from the nuclear region, our data are not of sufficient spectral 
resolution to constrain more than one velocity component. 

We check the fluxes in Ha between our net Ha image and 
our VIMOS IFU in a region on the inner south-east filament. 
The region is of size 19.5x23 arcsec centred on 14^07 m 30 s .2, 
-27°01'04.8". Our fluxes are 5.1 xlO -14 erg cm -2 s -1 and 
5.7 x 10 -14 erg cm -2 s -1 in the narrow-band image and IFU im- 
age respectively. The fluxes agree within ~ 10 per cent. The error 
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Rest wavelength / A 

Figure 5. The optical spectrum of the BGG nucleus in Abell 3581. 




30.4 14:07:30.0 29.6 



Figure 6. The line-of-sight velocity of the narrow line component of the optical nebulosity (left) from a fit to the strong [N n]A6583 line and the FWHM 
velocity width (right) after correction for the instrumental resolution (98 km s _1 , FWHM). Colour bars are in units of km s _1 and the velocity zero point 
is taken to be 6450 km s _1 , the velocity of the Ho; nebulosity at the position of the nucleus (14^07 m 29 s .8, — 27°01'06.4", marked with an X). Only pixels 
where the [N n]A6583 emission line was detected above Act are shown. 



on the integrated flux from our ID L fit ~ 7 per cent and a 1 per 
cent continuum flux variation between the ON/OFF filter images 
translates to ^10 per cent net flux change for A3581. Our fluxes 
are therefore consistent. 

The left hand plot of Fig. 7 shows the ratio of the strong [N 
II]A6583 emission line to the narrow component of the Ha emis- 
sion. In the centre of the galaxy this is observed to peak at the nu- 
cleus and to drop radially with distance from the AGN. Farther out 
onto the filaments the ratio is approximately unity. On the 'inside' 
of the filaments, the side closest to the radio emission, a larger ratio 
of [N II] to Ha emission is observed. 

The [O lll]A5007 emission line can only be observed in the 
eastern quadrants of the data cube; we have no information about 
the [O m]A5007 emission in the nucleus (right hand panel of Fig. 
7). The [O III]A5007 emission is over three times stronger in the 
easterly extension to the north of the nucleus, where the velocity 
widths are very high, compared to the southern-most regions of the 
easterly filament. A gradient of [O III] A5007 / Ha emission, similar 
to that observed in the [N ll]A6583 / Ha ratio discussed above, is 
observed across the more southern filament. The collisional line 
emission is enhanced relative to the recombination line emission 
on the 'inside' of the filament and close to the nucleus. 

Unusually strong emission lines of [N l]A5199, characteris- 
tic of the filament emission in BCGs, are detected in the extended 



filaments but no evidence is found for such emission in the nu- 
cleus (see Table. 2), we also detect He I emission in the extended 
filaments. We do not observe any evidence for coronal emission 
lines (emission lines with excitation temperatures of 10 5 — 10 6 de- 
gree gas) in our IFS data. Table. 2 gives emission line ratios com- 
pared to Ha for the three regions shown by crossed on Fig. 7. These 
are the nucleus, the north-easterly extension and farther out on the 
south-easterly extension. It is quite apparent that these ratios are 
different due to the comparative importance of various excitation 
mechanisms in these regions; this highlights the complicated na- 
ture of distinguishing excitation mechanisms in the central regions 
of BGGs/BCGs. 



5.1.4 Star formation 

A young star cluster is observed, projected on the dust lane, 2.5 
arcsec north of the nucleus (see Fig. 3). No other young clusters 
are observed in the filament system. Using our relative reddening 
map, and assuming the intrinsic dust is an obscuring screen, we cor- 
rect the observed fluxes in the broadband filters and determine the 
emission from the young star cluster. Using a GALEV evolutionary 
synthesis model (Kotulla et al. 2009) for a simple stellar popula- 
tion the F300X-F555W and F300X-F814W colours indicate an age 
<10 8 yrsandamass-10 4 M . 
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Figure 7. Left: [N ii]A6583/Hai ratio map of the filaments surrounding the two inner bubbles. The ratio is highest in the centre of the galaxy and smoothly 
drops at larger radius. A spur of higher [N n] is also seen north west of the centre where the velocity dispersion in the gas is larger. Right: The [O iii]A5007/Hai 
ratio map on the same spatial and intensity scale as the right hand image. The ratios are corrected for Galactic extinction but no correction for internal extinction 
has been made. The three crosses indicate regions (the nuleus, NE and SE) where spectra were extracted and the emission line fluxes detailed in Table 2. 



Table 2. Emission line fluxes (upper panel) and ratios (lower panel) in the 3 
regions indicated by a cross on Fig. 7. Each region was 2 pixels (~1.3 arc- 
sec) in radius. The units for the fluxes are x 10 — 15 erg cm - 2 s — 1 . 1 Errors 
in nucleus large - probably due to multiple components. 2 Spectral cover- 
age of VIMOS quadrants is uneven so [O m]A5007 emission can only be 
probed to the east of the nucleus. 3 Broader line than Hen due to not being 
able to resolve the doublet emission. 



Wavlength 


Nuc 1 


SE 


NE 


redshift 


0.021695±6E-6 


0.021746±lE-6 


0.02032±3E-5 


velocity width 


301±6 


209±1 


355±36 


[O m]A5007 


N/A 2 


0.4±0.02 


0.4±0.02 


[N i]A5199 




0.3±0.05 3 


0.09±0.02 3 


He I 5875 




0.1±0.02 




[O i]A6300 


1.7±0.4 


0.7±0.01 


0.2±0.03 


[O i]A6363 




0.2±0.02 




[N n]A6548 


2.7±0.1 


1.5±0.04 


0.3±0.01 


Ha 6563 


4.6±0.1 


3.1±0.01 


0.6±0.03 


[N n]A6583 


11.4±0.07 


4.5±0.02 


1.1±0.03 


[S n]A6717 


3.2±0.2 


1.5±0.02 


0.6±0.02 


[S n]A6731 


3.0±0.2 


1.2±0.03 


0.4±0.1 


Line / Ha. 


[O m]A5007 


N/A 2 


0.13 


0.67 


[N I]A5199 




0.10 


0.15 


He I 5875 




0.03 




[O i]A6300 


0.37 


0.23 


0.33 


[O i]A6363 




0.06 




[N n]A6548 


0.59 


0.48 


0.5 


Ha 6563 


1.0 


1.0 


1.0 


[N n]A6583 


2.48 


1.45 


1.83 


[S n]A6717 


0.70 


0.48 


1.0 


[S n]A6731 


0.65 


0.39 


0.67 



Abell 3581 has not been observed by Spitzer. However, the 
central galaxy is detected in all four WISE bands. The WISE All- 
Sky Source Catalog, measures an extrapolated flux at 24 /xm as 
8.0 mJy. Using equation 10 of Rieke et al. (2009), and assuming a 
Kroupa (2002) IMF, indicates the star formation rate (SFR) in the 
central galaxy is ~ 0.21 M yr _1 . The near- and far-ultraviolet 
emission from Galex indicates 0.3 M yr _1 and 0.1 M yr _1 
(these values assume no intrinsic reddening correction). 

The net Ha image of IC 4374 is shown in Fig. 1 . The total Ha 



flux is 6.9xl0~ 14 erg cm" 2 s~\ assuming [N ll]6583/Ha = 1.5 
and [N ll]6548/[NII]6583 = 1/3 (an average from our VIMOS IFS 
data - we note this does not cover the entire extent of the filament 
system). If the nucleus is excluded (a circle with a radius of 2.5 
arcsec), the total flux is 4.8 x 10 -14 erg cm -2 s _1 . Using the Ken- 
nicutt (1998) relation, SFR (Moyr -1 ) = 7.9xl0 42 (L H <* / ergs 
s _1 ), the SFR is ~ 0.4 M yr _1 , assuming no intrinsic extinction 
for the Ha emission. This is similar to the mass deposition rate in 
the cool X-ray gas (0.4 M yr _1 , Sanders et al. 2010), however, 
no UV emission is observed in the filaments. 



5.2 X-ray 

5.2.1 Morphology 

The X-ray surface brightness, temperature and metallicity maps are 
shown in the core in Fig. 8 and on a larger scale in Fig. 9. New 
outer structures are indicated in Fig. 10. As previously noted by 
Johnstone et al. (2005), who examined a 7165 s Chandra ACIS-S 
observation of Abell 3581, an X-ray emitting point source lies in 
the nucleus and two inner cavities exist; one to the east and one to 
the west of the nucleus. These inner cavities coincide with relativis- 
tic synchrotron emitting gas observed as radio bubbles in 1.4 GHz 
radio emission. The power in these inner 1.4 GHz radio bubbles is 
2.5 x 10 23 W Hz" 1 (Johnstone et al. 2005). 

Fig. 8 reveals that the low volume filling, soft X-ray gas, de- 
tected by Sanders et al. (2010), lies in filaments. There is excellent 
spatial correspondence between the coolest, highest density X-ray 
emitting gas and the Ha emitting filaments. Clear structures are 
seen in the temperature and metallicity maps suggesting little mix- 
ing is occurring in the centre of Abell 3581. 

Additional structure is observed outside of the filaments; sur- 
face brightness edges to the east and south indicate two bays where 
a dearth of X-ray emission could be a signature of outer bubbles 
(see Fig. 10, white ellipses). An additional X-ray cavity (green el- 
lipse), detached from the nucleus and not observed before, is seen 
^45 arcsec to the west of the nucleus. We do not observe any radio 
emission, in our frequency bands, coincident with the depression 
in the X-ray image. The cool X-ray gas is observed to extend up to 
the new outer bubble but not farther. A clear edge is seen in both 
maps at the inner edge of the bubble. 

The large scale structure in the X-ray image shows a spiral 
surface brightness discontinuity (red curve) which wraps around 
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0.99 1.3 1.7 2 2.4 0.34 0.79 1.2 1.7 2.1 9.04e-11 1.35e-09 1.14^08 



Figure 8. From left: X-ray surface brightness image, X-ray temperature map, X-ray metallicity map and 0.5 keV emission measure map. The inner two 
depressions corresponding to the regions in which 1 .4 GHz radio emission is seen are circled on the X-ray surface brightness image. Contours of Hen emission 
are overlaid on the temperature, metallicity and emission measure maps. The colour bars are in units of keV, Z© and cm -5 arcsec -2 respectively. 




1.1 1.2 1.4 1.5 1.7 1.8 2 2.1 2.3 0-27 0.37 0.47 0.58 0.68 0.78 0.8S 0.99 1.1 



Figure 9. From left: X-ray surface brightness image, X-ray temperature map and X-ray metallicity map. Two depressions are observed which coincide with 
the 1 .4 GHz radio emission in the nucleus, an outer cavity is observed towards the south west beyond a sharp temperature and metallicity edge and a further 
surface brightness discontinuity, a likely cold front, is observed beyond this outer bubble. There is an apparent metallicity drop toward the centre, as seen in 
other similar objects, this could indicate the need for additional temperature components to the model (see discussion in section 5.2.2). 




Figure 10. The outer features. Left: Unsharpmasked Chandra X-ray image of Abell 3581. An image, smoothed with a gaussian with a =25 pixels has been 
subtracted from the original. The outer cavity is highlighted in green and the cold front in red. The two 'bays' to the east and south are indicated with dashed 
white ellipses. Middle: The X-ray surface brightness image after a smooth elliptical isophotal model was subtracted from the original data. The intensity scale 
was chosen to highlight the outer cavity and spiral. Right: The same regions overlaid on the X-ray temperature map. 
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5.2.2 X-ray metallicity 
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Figure 11. Deprojected profiles showing the temperature, density, pressure 
and entropy of a sector centred on the cool core and intersecting the spiral. 
The entropy and temperature of the gas in the spiral is lower than that of the 
gas between the spiral and the core. The black regions contain 6000 counts 
and are shown on the above plot with the outer bubble, spiral cold front and 
approximate width of the spiral also shown. The red regions contain 3500 
counts such that the statistics remain good while allowing two regions on 
the spiral. 



the nucleus in a clockwise direction. The X-ray gas is cooler, denser 
and higher metallicity within this structure and the entropy is lower 
than the surrounding gas, similar to the gas in the X-ray cool core 
and may be gas stripped from there by sloshing motions in the core. 
Fig. 1 1 shows the X-ray properties through a sector centred on the 
nucleus and intersecting part of the spiral structure. 



The X-ray metallicity maps (Fig. 8 and Fig. 9) show a high abun- 
dance spiral structure but a low abundance core to the galaxy. Most 
cool-core galaxy clusters and groups tend to show a central abun- 
dance peak (e.g. De Grandi & Molendi 2001). However, an offset 
peak and metal enriched spiral is also observed in the Centaurus 
cluster (Sanders & Fabian 2002), Ophiuchus Cluster (Million et al. 
2010) and HCG 62 (e.g. Gu et al. 2007; Gitti et al. 2010; Rafferty 
et al. 2012). HCG 62 has a central iron abundance of ^0.3 Z© 
which rises to solar at 20-30 arcsec (5.6 — 8.4 kpc) before dropping 
father out (Rafferty et al. 2012) while Centaurus has a central iron 
abundance of ^0.4 Z© rising to ~1.3 Z© 15 kpc from the nucleus 
(Sanders & Fabian 2002). Similarly, with a single temperature fit 
to the data we find an iron abundance of ^0.3 Z© in the central 
regions of Abell 3581 which rises to 0.8-1.0 Z© at a distance of 
^30 — 60 arcsec (13.5 — 26 kpc) coinciding with the spiral struc- 
ture. The metallicity then drops farther out into the ICM. 

Resonance scattering can effect the measured abundances in 
clusters and groups; however, the effect has not been able to fully 
explain the metallicity drops (e.g. Mathews et al. 2001; Gastaldello 
& Molendi 2004; Sanders & Fabian 2006). A discussion on the 
effect of LMXBs on the derived spectral properties is given in 
David et al. (2009, 201 1). The gas in groups and clusters, viewed in 
projection, is certainly multi-phase and another possibility for the 
metallicity drop in the central regions of the galaxy may be that a 
single temperature model is inadequate in representing the data. 

To examine whether the metallicity drop is due to the mis- 
representation of a multi-phase gas, we fit two temperature VAPEC 
models within Xspec (version 12.8.0a) to the X-ray emission in 
four regions within the cluster. These regions are shown in Fig. 12 
and the results of the fit are given in Table 3. The minimum counts 
in our regions are 5400 counts with 2 of the regions having greater 
than 12000 counts. The results of Table 3 show that the addition of 
a second temperature component in regions 1 and 2 significantly al- 
ters the abundances removing any evidence for a metallicity drop in 
the central regions. The spiral (region 3) has a similar best fit metal- 
licity as the core while beyond the spiral (region 4) the metallicity 
declines. 



5.2.3 Bubble energetics 

To determine the energetics of the new outer cavity we fit an el- 
liptical isophotal model to the cluster and subtract the median sur- 
face brightness at each radius. We then identify, by examining cuts 
through the resulting surface brightness image, an ellipse represent- 
ing the cavity. 

We define a sector of the cluster centred on the nucleus and 
extending between the edges of the ellipse (see Fig. 13) up to a ra- 
dius of 100 arcsec (44 kpc). We create annular regions within this 
sector. The width of each annulus is determined to be the width nec- 
essary to enclose at least 3000 counts; we do this to ensure a good 
fit to the spectrum. The inner 3 pixels are not used to avoid includ- 
ing the point source. We then fit both the projected and deprojected 
spectra within Xspec (as described in Section 3). 

Estimating errors on energetics from cavities is fraught with 
difficulties due to inferring a 3 dimensional shape for the bubble 
from a projected 2 dimensional image, not knowing the inclination 
of the rise direction to the observer, the 'fuzzy' edges of the bubble 
decrement and due to the shape of the overlying surface bright- 
ness profile of the cluster. The errors estimated and quoted in the 
following sections use the technique of O' Sullivan et al. (2011). 
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Table 3. Results from the multi-temperature fits to regions 1-4 in Fig. 12. 
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Figure 12. Regions 1-4 where the X-ray spectra has been extracted and a 
multi-temperature gas is fit to the data. The point source in region 1 has 
been removed. The results of the single temperature and two temperature 
fits are given in Table 3. 

Using this technique the line-of-sight depth is equal to the semi- 
minor axis and the uncertainty on the volume is given by assuming 
a minimum depth of half this value and a maximum equal to the 
semi-major axis. However, for the outer cavity, it must be noted 
that the length of the bubble in a direction perpendicular to the rise 
direction is given by the semi-major, not semi-minor axis of the 
ellipse; the errors quoted are unlikely to be reliable. 

The temperature, density, pressure and entropy profiles in the 
gas in this sector are shown in Fig. 13, with entropy given by 
K — ksT/nV 3 . At the position of the X-ray cavity (indicated by a 
green solid line) we detect a clear jump in the deprojected tempera- 
ture and a corresponding drop in the density. There is no significant 
change in the pressure which could indicate that the jump in tem- 
perature is due to the gas in the region around the bubble and not 



the outer layers of the cluster seen in projection with the cavity. 
However, the errors on the pressure are not constraining. 
We estimate the energy in the outer cavity using, 

E bubble = — ^PV, (1) 
7-1 

where 7 = 4/3 for a relativistic fluid, V is the volume of the cav- 
ity and P is the thermal pressure of the ICM, as described in Birzan 
et al. (2004) and Dunn et al. (2005). The volume of the bubble is de- 
rived making the implicit assumption that the bubble is a spheroid, 
two of whose axes are in the plane of the sky, with a plane of 
symmetry perpendicular to the rise direction. The ellipsoid volume 
V = (4/3)7i\Rii?|, where ifr = 3.8 kpc (8.8") is the length of the 
projected ellipsoid axis in the direction of the nucleus, in this case 
the semi-minor axis of the ellipse, and R2 = 8.4 kpc (19.3") is the 
length of the ellipsoid in a perpendicular direction, here given by 
the semi-major axis of the ellipse. Using this formalism, Ebubble — 
5.5(2.7 — 11.8 x 10 56 ) ergs. A similar treatment gives the energy 
in the inner cavities as Ebubble = 2.4(1.2 — 2.9 x 10 56 ) ergs, as- 
suming Rt = 3.6 kpc (8.4") and R 2 = 3.0 kpc (6.9"). Diehl et al. 
(2008) estimate similar energies for these inner bubbles with values 
of 3.49 x 10 56 and 3.37 x 10 56 for the inner bubble energies. 

The outer bubble rise time as measured by the sound speed, 
buoyancy and refilling timescales (Blanton et al. 2001; McNamara 
et al. 2000; Churazov et al. 2000) are 3, 4, and 5 x 10 7 years respec- 
tively and for the inner bubbles are and 0.7, 2 and 6xl0 7 years. 
Assuming the sound speed timescale for the inner bubbles and 
the buoyancy timescale for the outer bubble the cavity powers are 
Px-cav — 4.4 x 10 41 ergs s _1 for the outer cavity and P x -cav — 
1.1 x 10 42 ergs s _1 for the inner cavities. A simple scaling from 
the power in the 1.4 GHz radio emission (2.5 x 10 23 W Hz -1 ), us- 
ing the relation derived by Cavagnolo et al. (2010), gives P cav — 
3.7 x 10 43 ergs s -1 for the inner cavities. 

The point source removed X-ray luminosity, between 0.5- 
7 keV, in a circle of radius 9.4 kpc (21.5") which encompasses 
the inner bubbles, a circle of radius 28.3 kpc (65.0") which en- 
compasses the outer bubble and the 'cooling radius' (assumed to 
be where t coo i < 7.7 Gyr, r coo i ~ 57 kpc, Sanders et al. 2010), 
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Figure 13. Temperature, density and pressure profiles in a sector chosen to 
include the cavity and cold front. The widths of the annuli are determined 
by requiring a minimum of 3000 counts per annuli and are shown in the 
image above. The green solid line on the figure indicates the radius of the 
center of the outer cavity while the blue dashed line indicates the radius of 
the cold front. 



corrected for Galactic absorption are 2.1 xlO 42 , 7.8 xlO 42 , and 
1.1 xlO 43 ergs -1 respectively. 



5.3 Radio 

5.3.1 Morphology 

The morphology of the radio emission is shown in Fig. 14. The high 
frequency 1.4 GHz emission fills the inner cavities; these are sur- 



rounded by the ionised gas nebulae which coincides with the lowest 
temperature X-ray emission. This is consistent with the inflation of 
the bubbles leading to the displacement of the cool and warm ISM 
in the galaxy, due to the pressure of the relativistic particles. 

The left hand panel of Fig. 14 compares the morphology of 
the lower frequency radio emission to that of the high frequency 
radio emission and the X-ray gas. At lower frequencies the orienta- 
tion of the radio extension is rotated clockwise with respect to the 
high frequency emission. The rotation is in the same sense as the 
Ha filament (see Fig. 14, right hand panel) and the 614 MHz con- 
tours (green contours) show a spur of emission which follows the 
filament in the north-easterly direction. 



6 DISCUSSION 

A distinguishing feature of many BCGs in cool core clusters, where 
heating is required to prevent catastrophic cooling of the hot gas, 
is that they harbour intricate extended systems of low ionisation 
optically emitting gas. A significant body of work has studied the 
origin of these filaments and their relationship with star formation, 
AGN feedback and gas in other phases. The BGG IC 4374 in Abell 
3581 is one such system where it is clear that the structures in the 
hot, warm and cold gas are intimately related and are connected 
with the presence of the radio lobes. 

IC 4374 exhibits prominent dust lanes running north and south 
through the galaxy nucleus. Both these dust lanes have extensions 
to the east which trace the ionised gas filaments (Fig. 3). A young 
star cluster is observed in the northern wedge shaped dust lane but 
no other evidence for star formation, notably none in the filament 
system, is observed. 

Scaling from our E(B — V) map (Fig. 4) of relative intrinsic 
dust extinction and assuming a hydrogen column density of Nh = 
5.8 x 10 21 x E(B - V) and the Galactic gas to dust ratio, ~100, 
allows us to estimate the dust mass in the galaxy as 10 6 M with 
~10 5 M in the south-east filament. 

The morphology of the Ha+[N II] emission is matched very 
well by the coolest and highest density regions of the X-ray gas 
(Fig. 8). The optical filaments engulf the inner bubbles and are co- 
incident with the easterly extending dust lanes. The smooth velocity 
gradients of the ionised gas suggest the filaments extending around 
the bubbles are continuous. 

The south-west filament leads to an outer bubble. Cool, low 
metallicity X-ray gas also extends from the centre of the cluster in 
the direction of this bubble and there is a discontinuity in the prop- 
erties of the X-ray gas at a radius corresponding to the inner edge 
of the bubble. The lower frequency radio emission is rotated clock- 
wise with respect to the high frequency emission and is aligned 
along the axis of the outer, not the inner, bubble (see Fig. 14). 

A spiral of high metallicity, low temperature, low entropy X- 
ray gas, with a cold front at its outer edge, is observed to trail from 
the cool core. In the following sections we draw on these observa- 
tions to attempt to answer questions pertinent to the nature of the 
heating and cooling cycle in this system. 



6.1 Origin of filaments and dust lanes 

The filaments in Abell 3581 exhibit point symmetry but are not 
radial, varying direction as they extend out of the galaxy. We cannot 
comment on the extent of the directional change as this depends 
on the projection angle. The velocity structure of the filaments is 
smooth, consistent with the cool and cold gas being pushed aside 
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Figure 14. Left: X-ray temperature map overlaid with contours of radio emission. Red contours are from 244 MHz GMRT emission, green contours from 
614 MHz GMRT emission and the central blue contours are from L-band EVLA radio emission. Low temperature gas (darker colours, see Fig. 8 and 9) is 
observed at the centre of the system and stretches east and west from the nucleus towards the bubbles. Right: The same radio contours overlaid on the SOAR 
net Hen image. The beam sizes are shown in the top right hand corner of the image. 



due to the expansion of the inner radio bubbles and, importantly, 
both the eastern filament extensions (the blueshifted and redshifted 
sides of the filament) are accompanied by thin threads of dust. 

Initial dust grain formation requires high densities and low 
enough temperatures to permit the condensation of silicates and 
graphites; it therefore cannot form straight out of a dust-free ICM. 
However, Draine (2003) makes the point that our current knowl- 
edge of grain destruction time-scales and observed element deple- 
tions in the gas-phase highlight that much of the subsequent dust 
growth must happen through interstellar processing. Dust and poly- 
cylic aromatic hydrocarbons (PAH) have both been found in the ex- 
tended ionised gas filaments of other BCGs (e.g. Sparks et al. 1989; 
Johnstone et al. 2007; Donahue et al. 2011). 

The south-west filament ends at the discontinuity in surface 
brightness, temperature and metallicity delineating the inner edge 
of the outer bubble. The correlation observed between the cool gas 
and dust filaments, X-ray cavities and the radio emission and the 
symmetry of the two filament extensions is hard to reconcile within 
the framework of filaments cooling from the hot ICM in situ or 
being the result of a galaxy merger. 

It is therefore tempting to advocate that this outer bubble and 
another in the opposite direction are responsible for the displace- 
ment of cool gas seeding the extended filaments and an interaction 
between the cool uplifted material and the hot ICM is responsible 
for the cooler X-ray gas (e.g. see discussion in Werner et al. 2012) 
However, we note that by eye, we do not observe a second outer 
bubble. The X-ray surface brightness is not symmetrical about the 
nucleus; the surface brightness drops more sharply in the north-east 
of Abell 3581 than in the south-west (Fig. 15) which might inhibit 
our view of a second outer bubble. 

Crudely, we test the hypothesis that the corresponding north- 
eastern outer bubble is not observed due to the low surface bright- 
ness by examining the significance of a similar outer bubble (same 
size/shape), reflected through the nucleus of Abell 3581. 

To quantify the significance of a bubble we compare the 
counts in the 'bubble region' to regions on either side of the bubble 



and with the same shape and radius. For our detected outer bub- 
ble we find a ~10 per cent decrement in counts (^3010 counts 
compared with ^3390 counts in the undisturbed regions) corre- 
sponding to a ^6 a deficit. If we assume a similar size and shape 
bubble with point symmetry through the nucleus, we find a ^5 per 
cent decrement (^2119 and ^2196 counts) corresponding to ~1.5 
sigma. A 10 per cent deficit in counts in this low surface brightness 
region, similar to the deficit seen in the outer bubble, would still 
only amount to a 3 sigma detection of a broad bubble and is proba- 
bly not easily observed by eye. We may not expect to see a bubble 
with the same size and shape and at the same radius. The surface 
brightness drop indicates a lower gas density. We might expect a 
bubble to be larger in size in a less dense medium further decreas- 
ing the likelihood of observing it. We cannot rule out the possibility 
of another outer bubble being responsible for the gas uplift of the 
filaments towards the north-east. We also note that a depression in 
the surface brightness profile is observed at ^40 arcsec from the 
nucleus towards the north-east direction (see Fig. 15) which lends 
supporting evidence for another outer bubble. 

The emission lines in the central regions of the galaxy have 
very large measured velocity dispersions particularly towards the 
eastern inner cavity. Our data are unable to constrain two or more 
velocity components to the line emission but it is likely that mul- 
tiple velocity components contribute to the broad velocity features 
observed. Farage et al. (2012) suggest the broad blueshifted feature 
is indicative of an inflow, perhaps responsible for feeding the AGN. 
We suggest the line-of-sight velocities, blueshifted in one direction 
and redshifted in the other, are consistent with the expansion of a 
radio bubble pushing aside the cool gas as it expands. This cool 
gas may have existed in an extended structure due to outflows from 
previous bubbles. The action of the current inner radio bubbles is 
to sweep aside the cool gas and to expand in the direction of least 
resistance. The coincidence of both the north-east and south-east 
filaments with filamentary dust lanes supports the idea that the gas 
is being lifted/pushed away from the centre of the galaxy. 

High spatial resolution observations of the cold gas, for exam- 
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Figure 15. X-ray surface brightness profile towards the outer X-ray cavity 
(south-west) and in a direction reflected through the nucleus (north-east). 
The X-ray emission is not radially symmetric; in the north-east direction the 
cluster surface brightness drops much more sharply than towards the south- 
west. The outer cavity in the south-west is at a distance of ~55 arcsec. 

pie with ALMA or SINFONI IFU, would shed light on the presence 
of an inflow or outflow in the nucleus. 



6.1.1 Origin of the gas central to the galaxy 

Having argued that the gas in the filaments originates in gas uplifted 
from the central galaxy the question then becomes, did the cool 
and cold gas in the galaxy originate from the cooling hot ICM or 
through stellar mass loss? This is a complicated question and the 
answer is probably both but whether a dominant mechanism is at 
play is not known. 

Abell 3581 has a total Ha flux of 6.9x 10" 14 erg cm -2 s _1 ; 
this translates to a mass, M H + = 2.6 x 1O 6 M assuming n e ~ 
100 cm -3 and that the effective recombination coefficient of Ha 
isa eff = 1.17 x 10" 13 cm 3 s _1 (Osterbrock & Ferland 2006). 
CO (2 — 1) emission has also been detected in the central galaxy 
(Edge, private communication). We assume ~10 8 M , of molecu- 
lar gas and estimate the observable dust mass from our E(B — V) 
map as 10 6 M . 

Sanders et al. (2010) find a mass deposition rate of about 0.4 
solar masses a year in the central regions of Abell 3581 and a cool- 
ing time for this gas of ~ 10 8 yr. The coolest X-ray gas is mor- 
phologically similar to the optical filaments. If these filaments are 
seeded by being dredged up out of the galaxy by the bubbles then 
the outer filaments are approximately the age of the outer bubble, 
~ 5x 10 7 years. In this time, making the extreme assumptions that 
100 per cent of the gas cools and that the mass deposition rate has 
not changed, the X-ray contribution to the filaments could be up 
to 2xl0 7 M , easily accounting for the mass of cool gas in the 
filaments but an order of magnitude shy of the expected cold gas 
content of the central galaxy. If the filaments are currently growing 
in mass, some dusty gas must have existed to seed the filaments 
in the first place. It should be noted that we don't know the exact 
morphology of the cold molecular gas; however, cold molecular 
gas has been found in many extended filament systems (e.g. Hatch 
et al. 2005; Johnstone et al. 2007; Oonk et al. 2010; Mittal et al. 
201 1 ; Salome et al. 201 1 ; Donahue et al. 201 1 ; Werner et al. 2012). 

The SFRs given by WISE infra-red and by GALEX near- and 
far-ultraviolet emission are ~ 0.21 M yr -1 , 0.3 M yr -1 and 
0.1 M yr -1 respectively, close to the current mass deposition rate 
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Figure 16. The time evolution of mass loss for a simple stellar popula- 
tion (SSP) which has been born at time t = using the models of BC03 
(Bruzual & Chariot 2003). The fraction of the stellar mass returned is very 
initial mass function (IMF) dependent but only midly dependent on metal- 
licity. 



measured from the cool X-ray gas. However, the morphologies of 
these two components are very different. The star formation is cen- 
tred on the nucleus and a single bright young star cluster, and is 
roughly symmetrical, while the cool X-ray gas is asymmetrically 
distributed and much of it is associated with the filaments. 

Fig. 16 shows that in the lifetime of a BGG a significant frac- 
tion of its stellar mass will be returned to the ISM and its various 
gas phases. The evolved stellar population is currently contributing 
~1 M yr -1 (estimate for a 10 11 M galaxy). This is comparable 
to the cooling rate of the hot gas and the SFR, and as such the stel- 
lar mass loss may be an important contributor to the current cool 
and cold gas being generated in the galaxy. 

We are unable to conclude whether one mechanism dominates 
over the other; both stellar mass loss and the cooling gas may be 
providing similar cool gas quantities currently. The stellar mass 
loss in total, over the lifetime of the population, must be signifi- 
cant, however, how much of the mass is maintained in the cool/cold 
gas phases is not known. How X-ray mass deposition rates and the 
heating of the X-ray gas has varied with time is also not known. 
The more abundant cool and cold gas in BCGs within galaxy clus- 
ters with cool cores, as opposed to non-cool core galaxy clusters, 
suggests that the high pressure environment must play a role. 
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6.2 Filament excitation 

Clues about the filament excitation mechanisms can be gleaned 
through Fig. 7. Near the nucleus there is both a gradient in the [N 
II]A6583 / Ha emission line ratio, decreasing rapidly from the nu- 
cleus (dropping by a factor of 2 within 5 arcsec), and a gradient 
from the inside to the outside of the nearby filaments (dropping 
by a factor of ~1.5 across the filament width). Farther out in the 
filaments the ratio shows less variation. 

At a fixed metallicity and density the [N II] A6583 / Ha ratio is 
essentially an indicator of the amount of heating in the gas per hy- 
drogen ionization. However, collisional excitation may be more ef- 
ficient at higher densities, and increases in the nitrogen abundance 
in the gas will increase the [N ll]A6583 / Ha ratio. The observed 
gradient could be indicative of an extra heating mechanism, such as 
photo-ionisation by the AGN or a weak shock or interaction driven 
by the AGN outflow, or a metallicity or density gradient in the gas. 

As we suspect the gas originates in the galaxy and has been 
dragged out and pushed aside by the inflating radio bubbles, strong 
metallicity gradients over small regions of the filament are unlikely. 
[O m]A5007 is typically weak in the extended filaments in BCGs. 
Indeed it is too weak to observe in our outer filaments, but is strong 
in the filament sections close to the nucleus and where the veloc- 
ity widths of the lines are large. We observe a strong gradient (> 
factor 3 increase) in [O III]A5007 emission from the 'inside' of the 
southern filament, adjacent to the radio bubble, to the 'outside' and 
an accompanying strong gradient in the ratio of [O III] A5007 to Ha 
emission. Unfortunately the wavelength range of our spectrum does 
not include [O lll]A4363 emission so we are unable to estimate the 
kinetic temperature variations across the filament. 

[S ll]A6716 / [S ll]A6731 is a diagnostic of the electron den- 
sity. Ratios > 1.3 are typical of the low density limit (<100 cm -3 ) 
while the high density limit is characterised by ratios of < 0.5 
(>1000 cm -3 ). Typically, extended filaments have been found to 
lie in the low density limit (e.g. Hatch et al. 2006) and this is con- 
sistent with our observed values. 

We observe emission lines of [N l]A5199 in the filaments. This 
emission is usually weak in gas with a well defined kinetic temper- 
ature. This is due to the excitation potential being similar in mag- 
nitude to the ionisation potential of H°. There is often little N° 
present in warm gas and so little [N l]A5199 emission (see Ferland 
et al. 2009, 2012 for further discussion of [N l]A5199 excitation 
mechanisms). A mechanism which produces supra- thermal parti- 
cles in the gas and can create some level of ionisation even within 
cold gas, gas which is multi-phase on a microscopic level, would 
be successful at reproducing these line ratios (Ferland et al. 2009). 
Such a mechanism, provided by the surrounding hot ISM/ICM, has 
been suggested by Ferland et al. (2008, 2009); Fabian et al. (201 1). 
Additionally, Churazov et al. (2013) have suggested non-thermal 
particles created in situ in the magnetically supported filaments 
may contribute to the emission line spectra. 

The variation in the emission line ratios shown in Table 2 high- 
lights the complicated nature of the filament excitation mechanisms 
in regions which are not isolated, such as those close to the nucleus 
or close to inflating radio bubbles. Clearly several excitation mech- 
anisms are at play in the inner regions of this filament system and 
disentangling them is not trivial. It may not be informative to study 
the emission line ratios in large regions of the centres of these sys- 
tems; choosing specific 'clean' regions where one mechanism dom- 
inates is preferable. 



6.3 Jet precession/bulk motions of ICM 

The current bubble axis in IC 4347 has east-west point symme- 
try through the nucleus. However, the outer bubble appears twisted 
~20 degrees south from the westerly direction (in projection); 
a twist which is also reflected in the trailing Ha filaments and 
the high-to-low frequency radio emission. Lower frequency radio 
emission probes less energetic plasmas, therefore older bubbles. 
The effect of bubble wobbling is unlikely to account for the change 
in the orientation of the radio emission axis as this need not affect 
both the bubbles in the same way. 

Bubble rise directions do not appear to be consistent; some 
systems exhibit directional changes between outbursts but, for the 
most part, maintain point symmetry (e.g. NGC1275, Fabian 2012). 
Dunn et al. (2006) suggested a precessing jet may be responsi- 
ble for the directional change in the three most recent outbursts of 
NGC 1275. A motivating feature for this model is that the filaments 
are largely radial. The bubble direction may also have changed due 
to turbulence or gas motions within the nuclear region and subse- 
quently risen in the 'easiest' direction - that of the lowest pressure 
- however, it is not obvious that this would lead to bubbles with 
approximate point symmetry. If the precession model is an expla- 
nation for the Abell 3581 bubbles, assuming the bubbles are rising 
in the plane of the sky, the measured angle between the inner and 
outer bubbles in Abell 3581 gives a precession rate of 6 degrees per 
10 7 yrs. 

NGC 5044 exhibits bubbles which are rising in various direc- 
tions in the cluster (David et al. 2009, 2011). The authors suggest 
that the nearly isotropic distribution of the bubbles is a result of the 
group 'weather' from AGN driven turbulence generated by many 
small outbursts or from the sloshing of NGC 5044 with respect to 
the group potential (e.g. see simulations of Briiggen et al. 2005; 
Heinz et al. 2006). Finally, a beautiful example of linearly rising 
bubbles is shown by Randall et al. (2011) in their X-ray images of 
NGC 5813. 

Some clues as to the directional change in Abell 3581 are 
given by observations of the same twist in the optical and radio 
emission and the large-scale spiral structure observed in the hot X- 
ray gas beyond the outer bubble (Fig. 9). Fig. 13 shows evidence 
for a jump in temperature, a density dip and approximate pressure 
equlibrium coincident with the outer edge of the spiral, consistent 
with a cold front. 

Similar structures have been observed in many X-ray images 
of cluster cool cores (e.g. Churazov et al. 2000; Fabian et al. 2000; 
Mazzotta et al. 2001; Markevitch et al. 2001, 2003; Mazzotta et al. 
2003; Churazov et al. 2003; Dupke & White 2003; Clarke et al. 
2004; Sanders et al. 2005, 2009). In relaxed clusters, the current 
paradigm is that these cold fronts are caused by gas sloshing within 
the cluster cool core. The sloshing may be a result of a disturbance 
caused by the in-fall or fly-by of a sub-cluster or alternatively a 
disturbance associated with energy injection into the ICM by the 
AGN (for a review see Markevitch & Vikhlinin 2007). 

Jones et al. (2002), Baldi et al. (2009) and David et al. (2009) 
have found similar spiral features in NGC 4636 and NGC 5044 re- 
spectively with the distinguishing feature that these spirals appear 
to have two arms. Simulations of cluster sloshing have thus far pro- 
duced one-armed spirals. In NGC 5044 David et al. (2009) suggests 
the the group weather has been responsible for perturbing the buoy- 
ant bubbles rise paths causing them to drag cooler, metal rich ma- 
terial out of the galaxy in two spirals; one emanating from each jet. 
The 'weather' is perhaps still generated by bulk sloshing motions 
but the cooler gas is now dredged up from the nucleus rather than 
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cooling from the core. NGC 4636 differs in that its spiral arms are 
hotter, not colder than the surrounding medium suggesting these 
are the result of shocks driven into the ICM due to the inflation of 
the bubbles. 

In Abell 3581 the X-ray data appears to favour a single spiral 
structure. The slow rotation in the radio emission and the direc- 
tional change along the filaments indicates the rise direction dif- 
ferences are not due to jet precession and that some other process 
changed the bubbles direction after it had left the nucleus. This is 
unlikely to be simply the path of least pressure as there is a ro- 
tational symmetry to the directional change in the filaments. The 
projected twist direction is clockwise and is the same direction in 
which the spiral unravels. The stirring of the ICM from some pre- 
vious disturbance could be responsible for this. 

6.4 Cooling and heating 

Active AGN heating is ongoing in the core of Abell 3581. Bubbles 
are being produced and rising into the ICM dragging with them 
multi-phase filaments. The cool and cold gas and dust phases are 
presumably interacting with the hot surrounding medium resulting 
in the cool, low volume filling, X-ray emission. 

The X-ray luminosity in the core (at a radius of the inner bub- 
bles) is extremely well matched by the powers of the inner bubbles, 
L x ~2.1xl0 42 ergs" 1 and P cav ~ 2.2 x 10 42 ergs s" 1 re- 
spectively. However, assuming the outer bubble to the west of the 
nucleus is part of a pair, the outer and inner cavity powers, V cav ~ 
3.1 x 10 42 ergs s _1 , are more than a factor of two lower than the 
X-ray luminosity in a circular region with a radius of the projected 
distance of the outer cavity (Lx ~7.8x 10 42 erg s _1 ) and signifi- 
cantly smaller than the X-ray luminosity within the 'cooling radius' 
(assumed to be where t coo i < 7.7 Gyr, Lx ~1.1 x 10 43 erg s _1 ). 
If the bays observed to the south and east are additional bubbles, the 
total cavity heating within this region may be significantly larger. 

The high metallicity, low entropy spiral of gas eminating from 
the core is indicative of ongoing stripping of the X-ray cool core. 
This stripping is distributing low entropy material over a larger ra- 
dius in the cluster in a manner reminiscent of that observed in the 
much more massive Ophiuchus Cluster (Million et al. 2010). The 
mass of the cool core inside a radius of 30 arcsec (^13 kpc) is 
~10 10 M . Taking the density from Fig. 11 and a conservative 
cylindrical volume representative of the brightest part of the spiral, 
with dimensions 2 arcminutes in length by 10 arcseconds in radius, 
we estimate the mass of the stripped spiral arm to be ~8 x 10 8 M 
(making the assumption of 100 per cent filling factor) or ~8 per 
cent of the current cool core. Assuming the sound speed of the gas 
^400 km s -1 , the time taken for gas to travel the 2 arcmins from 
the core to the end of the cylindrical volume in which we have 
measured the gas mass is ~1.3 x 10 8 yr. Thus, stripping of this rate 
could in principle destroy the cool core in 1-2 Gyr, although ongo- 
ing cooling (the cooling time in this region is comaprable, ~10 9 yr) 
will also replenish it. 



7 CONCLUSIONS 

Abell 3581 has a rich multi-phase phenomenon not common 
in X-ray groups. There is evidence for multiple 'radio-mode' 
outbursts of activity from the central AGN and it is clear the 
resulting bubbles are interacting with the cold, cool and hot gas in 
the system. We summarise our main conclusions below: 



• The BGG, IC4373, exhibits bright cool ionised gas filaments 
extending over 30 arcsec (~13 kpc) from the nucleus. These 
filaments are coincident with dust and soft X-ray emission. The 
filaments engulf the inner east-west bubbles and then change 
direction extending towards older cavities in the south-west and 
north-east. 

• Based on the morphology of the X-ray cavities, radio emission 
and cool gas filaments and on the filament kinematics we suggest 
gas uplift is the origin of the ionised gas nebulosity. 

• This conclusion is supported by the presence of the dust in the 
filaments which is seeded by stellar mass loss in the centre of the 
galaxy. 

• The line emission close to the nucleus and coincident with the 
direction of the jet has very high velocities indicative of an outflow. 

• The BCG has a SFR of about 0.2 - 0.3 M yr -1 , as measured 
by WISE and GALEX. A young blue star cluster is also revealed 
in the HST data. The SFR is comparable to the X-ray mass 
deposition rate (0.4 M yr -1 ) from the XMM RGS spectral 
analysis (Sanders et al. 2010). 

• Farther out on the filament system, the emission line ratios are 
similar to those seen in the extended filaments of other BCGs 
exhibiting characteristically low [O lll]A5007 emission and 
relatively strong [N l]A5199 emission. It is probable that many 
excitation mechanisms are occurring in the 'messy' inner regions 
of the galaxy while the outer 'cleaner' regions are consistent with 
an excitation mechanism which produces some level of ionisation 
even in cold gas, such as collisions with supra-thermal particles 
from the hot ISM/ICM. 

• The presence of the cold front and the directional change of the 
bubbles and filaments indicates that the group weather, driven by 
sloshing, may be controlling the rise direction of the bubbles. 

• The spiral cold front exhibits higher metallicities and lower 
entropy than the surrounding gas and is similar in its X-ray 
properties to the cool core. The spiral could be material, stripped 
from the core due to sloshing motions. 
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